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Abstract 


The  “net  shape”  HIP  approach  is  based  on  computer  programs  with  embedded 
engineering  models  of  powder  consolidation,  shrinkage  and  surface  formation  that  use  the 
computer  aided  design  (CAD)  of  the  part  needed  and  of  the  HIP  tooling.  Densification  pattern 
during  HIP,  which  is  the  “key”  to  “net  shape”  technology  is  determined  by  the  capsule  and 
powder  materials  rheology  and  is,  at  the  initial  stages,  controlled  by  the  capsule’s  plastic 
stiffness. 

However,  dilatometric  experiments  on  HIP  of  capsules  with  powder  revealed  that  even  at 
the  final  stages  of  densification  when  the  capsule  is  pliable,  deformations  are  not  uniform  and 
follow  the  dominating  radial  or  axial  component  though  the  pressure  is  isostatic  and  of  uniform 
density.  Deformation  pattern  is  determined  by  a  form  of  deformation  history,  such  that  some 
“anisotropy”  generated  during  this  history  becomes  inherent  and  this,  in  turn,  controls  the 
process  of  shrinkage  even  under  uniform  pressure  of  HIP. 

Following  this  approach,  the  equations  describing  the  material  model  during  HIP  have  been 
modified  and  appropriate  experiments  carried  out  for  their  parametrical  identification. 

Anisotropic  behavior  of  powder  material  during  “net  shape”  HIP  consolidation, 
discovered  during  special  experiments  with  porous  samples,  was  demonstrated  experimentally 
and  by  introducing  an  anisotropy  module  within  the  entire  HIP  process  numerical  model. 

Basic  Hypotheses  and  Equations 

The  advancement  towards  “net  shape”  HIP  is  provided  by  the  process  modeling  and  HIP 
tooling  design  .  The  proposed  steps  to  highly  desirable  “net  shape”  products  are  improvements  to 
existing  process  models  used  to  design  HIP  tooling.  The  more  advanced  modeling  must  fully 
account  for  the  evolution  of  material  properties  during  HIP  cycles.  In  other  words,  the  powder 
rheology  may  depend  on  the  inherent  anisotropy  caused  by  the  evolution  of  pores  during 
consolidation. 

These  mechanisms  provide  additional  control  of  capsule/powder  shrinkage  during  HIP. 

To  introduce  inherent  anisotropy  caused  in  HIP  triggered  by  initial  difference  in 
capsule  stiffness  and  then  developed  due  to  the  pore  evolution,  it  is  necessary  to  accept  the 
following  basic  hypotheses. 


1.  Hypothesis  Nl. 


where  £tj 


There  exists  a  plastic  potential ,  i.e.,  there  exists  a  yield  surface, 


and  fik  -are  some  kinetic  parameters. 


/(<V/**)  =  °, 


2.  Hypothesis  N2. 

“Distribution  A.  Approved  for  public  release;  distribution  unlimited.” 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

APR  2005 


2.  REPORT  TYPE 


3.  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

HIP  Modeling  Methodology  Based  on  the  Inherent  Process  Anisotrophy 


6.  AUTHOR(S) 

Victor  Samarov;  Vassily  Goloveshkin;  Charles  Barre 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Synertech  P/M,  Inc., 11652  Monarch  St., Garden  Grove, CA, 92841-1817 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5a.  CONTRACT  NUMBER 
5b.  GRANT  NUMBER 
5c.  PROGRAM  ELEMENT  NUMBER 
5d.  PROIECT  NUMBER 

OSDB 

5e.  TASK  NUMBER 

R4PT 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

The  "net  shape"  HIP  approach  is  based  on  computer  programs  with  embedded  engineering  models  of 
powder  consolidation,  shrinkage  and  surface  formation  that  use  the  computer  aided  design  (CAD)  of  the 
part  needed  and  of  the  HIP  tooling.  Densification  pattern  during  HIP,  which  is  the  "key"  to  "net  shape" 
technology  is  determined  by  the  capsule  and  powder  materials  rheology  and  is,  at  the  initial  stages, 
controlled  by  the  capsule’s  plastic  stiffness.  However,  dilatometric  experiments  on  HIP  of  capsules  with 
powder  revealed  that  even  at  the  final  stages  of  densification  when  the  capsule  is  pliable,  deformations  are 
not  uniform  and  follow  the  dominating  radial  or  axial  component  though  the  pressure  is  isostatic  and  of 
uniform  density.  Deformation  pattern  is  determined  by  a  form  of  deformation  history,  such  that  some 
"anisotropy"  generated  during  this  history  becomes  inherent  and  this,  in  turn,  controls  the  process  of 
shrinkage  even  under  uniform  pressure  of  HIP.  Following  this  approach,  the  equations  describing  the 
material  model  during  HIP  have  been  modified  and  appropriate  experiments  carried  out  for  their 
parametrical  identification.  Anisotropic  behavior  of  powder  material  during  "net  shape"  HIP 
consolidation,  discovered  during  special  experiments  with  porous  samples,  was  demonstrated 
experimentally  and  by  introducing  an  anisotropy  module  within  the  entire  HIP  process  numerical  model. 


15.  SUBIECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

8 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


The  shape  (the  equation)  of  the  yield  surface  depends  only  on  the  current  values 
of  accumulated  defonnations,  not  on  the  process  history  (trajectory).  It  means  that  the  surface 
equation  is  as  follows:  f(<jiJ,ej)  =  0 . 


3.  Hypothesis  N3. 


The  initial  state  of  the  powder  material  is  isotropic.  It  means  that  in  the  principal 
system  axes  of  deformation  tensor  the  following  equation  is  valid:  /  =  /(c rjj,s1,s2,,s3)  =  0  . 

We  also  presume  that  the  equation  for  the  surface  can  be  presented  as  a  quadratic 

form. 


Assuming  all  of  this,  the  equation  for  the  yield  surface  can  be  presented  as: 

f  =  gl&l  +  §2Cr2  §3^3  +  £12^1^2  §13^1^3  S^23Cr2Cr3  + 

7  2  .  7  2  ,  7  2 

^12^12  +  ^13^13  +  ^23^23 

where  the  functions  g,  d  are  the  functions  of  s1,s2,s 3 . 

Let  us  consider  an  experiment  on  axisymmetric  upsetting  of  a  cylindrical 
specimen  with  the  density  p  along  the  axis  “3”. 

Then  s,  =  s2  =  s  ,  s2  =  s ,  2s  +  e  =  3e, 

Then  the  yield  surface  equation  can  be  presented  as  follows: 
g(s,p)al2  +  g(s,p)a22  +  g(s,  p)a32  +2i//(s,  p)a1a2+2i//(s,  p)axa3  + 

+2  y/{s,  p)(J2 ct 3  +2d(s,  p)oX2  +2 d(s, p)cu2  +2d(s, p)<J12~  =  a  2 

Then,  if  we  consider  upsetting  a  specimen  cut-off  at  the  angle  cp  to  the  axis  (1)  in  the  plane 
perpendicular  to  the  axis  (3)  the  results  should  not  depend  on  the  value  of  (p . 

Let  us  designate  upsetting  stress  as  cr . 

Then:  cr3  =  0  , 


^31  *-^32  0  , 


cr,  =cr  cos  cp , 


a2  =<j  sin2  cp  , 
cr12  =  cr  cos  ^  sin  . 

If  we  substitute  these  values  in  the  equation  for  the  yield  surface,  we  shall  obtain  the  following: 
cr2 [g(£,/?)(cos4  cp  +  sin4  cp)  +  2y/{s, /?)sin2  ^cos2  cp  +  2d{s, /?)sin2  cpcos2  cp\=<j  2 
or 

o' 2  [g(s,  p)  +  2  sin ^ 2  cp  cos ^ 2  cp(y/{s,  p)  +  d(s,p)~  g(e,  /?))]  =  <js 2 


Due  to  the  axial  symmetry,  the  value  of  cr  does  not  depend  on  cp . 
Then:  d(e,  p)  =  g(s ,  p)  -  y/{s,  p) 


Where 


2s  +  s  =  3e, 
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Plan  of  Experiments 
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The  following  experiments  were  carried  out  to  detennine  and  characterize  the 
functions.  There  will  be  two  types  of  experiments:  experiments  with  the  “isotropic”  material  (#1 
and  #2)  and  experiments  with  “non- isotropic”  material  (#3,  #4  and  #5). 

1.  Experiment  #1. 

HIP  in  a  thin-walled  symmetrical  capsule  with  the  deformation  tensor  being  close 
to  the  spherical  one.  Then, 


As  crj  =  a2  =  a3  =  P(p) ,  the  plasticity  criterion  will  give  the  following  relation  for  the  two 
functions: 

^2 

g(e1,p)  +  2i//(el,p)  =  jp 


2.  Experiment  #2. 

Uniaxial  deformation  (upsetting)  of  the  samples  obtained  in  Experiment  #1. 
Then,  cq  0,  <r2  =  cr3  =  0  and  as  a  result,  we  get  the  second  relation: 

g(et,p)  =  ^ 

A 

and  can,  therefore,  detennine  the  values  for  both  functions  in  the  point  (p,ei ). 

The  value  for  the  third  function  in  this  point  ( p,e{ )  can  be  derived  if  we  consider 
upsetting  of  this  isotropic  material  along  the  axis  other  than  “1”  which  leads  to  the  following 
relation: 

d(ex,p)  =  g(el,p)-y/(el,p) 

3.  Experiment  #3. 


HIP  in  a  thick  walled  capsule  leading  to  practically  uni-directional  deformation 
along  the  axis  “1”.  For  this  deformation  model 


—  0,  —  f?2 


A 


1 


P 


Letting  the  axial  stress  be  oa ,  and  using  the  equation  for  the  yield  surface  and  the  condition  that 
both  of  the  deformations  in  the  directions  orthogonal  to  “1”  are  equal  to  zero,  we  can  get  the 
following  relation 


Po 

for  the  functions:  g  =  g(s,p) ,  y/  =  y/(£’P)  in  the  two  points:  £=0,s  =  e2  =  —  -1 


P 


g(e2,p)~ 


2^2(0  ,p) 


<j. 


g(0,p)  +  i//(e2,p)  cr: 


4.  Experiment  #4. 
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Uniaxial  deformation  (upsetting)  of  the  samples  obtained  in  Experiment  #3  in  the 
direction  of  the  maximal  deformation. 

If  the  flow  stress  of  upsetting  is  si ,  then  the  equation  for  the  yield  surface  will 
provide  the  second  relation,  enabling  us  to  determine  g  =  g(e2,p): 

g(e2,p)  =  — 

si 

The  ratio  between  the  defonnations  in  the  direction  “1”  and  orthogonal  direction  (an  analogue  of 

8 '22 

the  Poisson  coefficient)  =  —  will  give  the  third  relation  between  the  functions: 

£w 


U(0,A> 

g(e2,P) 


5.  Experiment  #5. 

Uniaxial  deformation  (upsetting)  of  the  samples  obtained  in  Experiment  #3  in  the 
direction  orthogonal  to  the  maximal  deformation. 

If  the  flow  stress  of  upsetting  is  s2 ,  it  will  give  the  fourth  relation  for  the  functions: 

g(0,P)  =  ^t 

S2 

As  a  result,  we  can  obtain  the  values  of  two  functions  in  two  more  points  in 
addition  to  the  values  detennined  after  Experiments  #1  and  #2. 

For  the  third  function,  due  to  isotropic  properties  in  the  directions  orthogonal  to 
the  maximal  deformation,  we  get  the  relations  in  the  points 
d  =  d{e2  ,p) 
and 

d{e2,p)  =  g{b,p)-y/{e2,p) 

As  a  result  of  the  treatment  of  experimental  data,  we  get  the  values  of  the  first  two 
functions  in  three  points  ( s  =  0,  s  =  ex ,  s  =  e2),  and  for  the  third  function,  in  two  points 
(s  =  ex ,  £  ~  e2 ),  for  a  given  value  of  density  reached  during  HIP,  i.e.,  to  reveal  the  effect  of 
deformation  “anisotropy”. 

Carrying  out  experiments  for  several  HIPed  densities,  it  becomes  possible  to  build 
approximations  for  these  three  functions  g  =  g(s,p) ,  \j/  =  y/(c,p) ,  d  =  d(t\p) . 

For  the  isotropic  HIP  model,  the  functions  of  the  plasticity  criterion  depend  on 

density  only. 

Actually  “anisotropy”  introduced  in  the  proposed  model  is  based  on  the 
assumption  of  the  existence  of  the  “areal  density”  and  is  determined  by  the  values  of  density  and 
of  the  normal  deformation  in  the  principal  axes. 


Development  of  the  Database  for  Rheological  Coefficients  of  the  Anisotropic  Model 

Five  special  sets  of  experiments  with  isostatic  and  anisotropic  loading  were  envisaged  to 
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the  values  of  the  functions  g  =  g(s,p) ,  y/  =  y/{s,p) ,  d  =  d(s,p)  and  of  the  influence 
coefficient  [5  for  several  HIPed  densities. 

Titanium  powder  alloy  Ti-6A1-4V  was  chosen  as  a  candidate  to  gain  experience  on 
producing  complex  shaped  components.  Keep  in  mind  that  Ti-6A1-4V  only  showed  most  of  the 
distortions  and  hard-to-predict  features  for  the  titanium  alloy  rather  than  for  Ni-base  alloys. 

To  provide  uni-axial  defonnation  during  HIP  and  create  initial  anisotropy  of  the 
deformation  pattern,  special  thick  walled  capsules  with  the  increased  radial  stiffness  from 
mild  1018  steel  were  designed  and  manufactured  to  carry  out  parametric  identification. 

1.  Experiment  #1. 

HIP  in  a  thin-walled  symmetrical  capsule  with  the  deformation  tensor  being  close 
to  the  spherical  one. 


As  o"j  =  0%  =  o~2  =  P(p) ,  the  plasticity  criterion  will  give  the  following 
relation  for  the  two  functions: 

<j~ 

g(el,p)  +  2y/(el,p)  =  jp 

2.  Experiment  #2. 

Uniaxial  deformation  (upsetting)  of  the  samples  obtained  in  Experiment  #1. 
Then,  ^  0,cr2  =  a3  =  0  and  as  a  result,  we  get  the  second  relation: 

,  ,  o?  ' 

g(el,p)  =  — 

and  can,  therefore,  detennine  the  values  for  both  functions  in  the  point  ( p,e i ). 

The  value  for  the  third  function  in  this  point  (p,e{)  can  be  derived  if  we  consider 
upsetting  of  this  isotropic  material  along  the  axis  other  than  “1”,  which  leads  to  the  following 
relation:  d(eltp)  =  g(ex,p)~  y/ie^p) . 


3.  Experiment  #3. 


HIP  in  a  thick-walled  capsule  leading  to  practically  uni-directional  defonnation 
along  the  axis  “1”. 

For  this  deformation  mode 


£2  —  St,  —  0,  £■[  — 


Po 


1 


p 


Let  the  axial  stress  be  <Ja ,  then,  using  the  equation  for  the  yield  surface  and  the  condition  that 
both  of  the  deformations  in  the  directions  orthogonal  to  “1”  are  equal  to  zero,  we  can  get  the 
following  relation  for  the  functions:  g  =  g(s,p) ,  y/  =  y/(s,p)  in  the  two  points: 


s  =  0,  s  =  e2 


Po 


1 


P 
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g(e2,P) 


2^2(0,A> 
g{0,p)  +  y/(e2,p) 


4.  Experiment  #4. 

Uniaxial  deformation  (upsetting)  of  the  samples  obtained  in  Experiment  #3  in 
the  direction  of  the  maximal  deformation. 

If  the  flow  stress  of  upsetting  is  Sj ,  then  the  equation  for  the  yield  surface  will 
provide  the  second  relation,  enabling  us  to  determine  g  =  g(e2,p)  : 

g(e2,p)  =  — 

si 

The  ratio  between  the  defonnations  in  the  direction  “1”  and  orthogonal  direction 

6*22 

(an  analogue  of  the  Poisson  coefficient)  /(  =  —  will  give  the  third  relation  between  the 

sn 

functions: 

_  y(0,p) 

1  g(e2,p) 


5.  Experiment  #5. 

Uniaxial  deformation  (upsetting)  of  the  samples  obtained  in  Experiment  #3  in 
the  direction  orthogonal  to  the  maximal  deformation. 

If  the  flow  stress  of  upsetting  is  s2 ,  it  will  give  the  fourth  relation  for  the 

functions: 

a1 

g(Q>p)  =  ^t 

u 

As  a  result,  we  can  obtain  the  values  of  two  functions  in  two  more  points  in 
addition  to  the  values  detennined  after  Experiments  #1  and  #2. 

For  the  third  function,  due  to  isotropic  properties  in  the  directions  orthogonal 
to  the  maximal  defonnation,  we  get  the  relations  in  the  points 
d  =  d(e2,p) 
and 

d(e2,p)  =  g(0,p)-i//(e2,p) 

As  a  result  of  the  treatment  of  experimental  data,  we  get  the  values  of  the  first  two 
functions  in  three  points  ( £  =  0,  c  =  e, ,  s  =  e2 ),  and  for  the  third  function,  in  two  points  (s  =  e j , 
e  =  e2 )  for  a  given  value  of  density  reached  during  HIP,  i.e.,  to  reveal  the  effect  of  deformation 
“anisotropy”. 

Carrying  out  experiments  for  several  HIPed  densities,  it  becomes  possible  to  build 
approximations  for  these  three  functions  g  =  g(s,p) ,  y/  =  y/(s,  p) ,  d  =  d(s,  p) . 

For  the  isotropic  HIP  model,  the  functions  of  the  plasticity  criterion  depend  on 

density  only. 
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Actually,  “anisotropy”  introduced  in  the  proposed  model  is  based  on  the 
assumption  of  the  existence  of  the  “areal  density”  and  is  determined  by  the  values  of  density  and 
of  the  normal  deformation  in  the  principal  axes. 

Figures  1  and  2  present  the  results  of  micro-structural  analysis  of  the  mounts 
taken  from  the  samples  oriented  in  different  orthogonal  directions.  These  photos  definitely  reveal 
the  existence  of  the  “areal  density”  phenomenon,  when  the  deformation  and  density  pattern 
provides  initial  anisotropy. 

The  results  of  the  compression  tests  were  processed  in  order  to  be  incorporated 
into  the  HIP  model  and  to  update  the  anisotropy  module  with  actual  data  received  from  the 
experiments. 


xlOO  xlOO 

Figure  1.  Typical  Micro-structures  of  a  Porous  Sample  in  the  Direction 
Orthogonal  to  Main  Deformation  Axis  (Average  Density  80%) 


xlOO  xlOO 

Figure  2.  Typical  Micro-structures  of  a  Porous  Sample  in  the  Direction  of  the 
Main  Deformation  Axis  (Average  Density  81%) 


CONCLUSIONS 

Anisotropic  behavior  of  Ti-6A1-4V  atomized  spherical  powder  during  HIP 
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Consolidation,  caused  by  the  initial  non-unifonn  capsule  stiffness,  was  revealed  during  special 
experiments.  The  values  of  the  flow  stress  ^  of  the  samples  cut  off  in  the  direction  of  the 
maximal  deformation,  processed  in  Experiment  #4,  are  about  30%  higher  than  the  flow  stress  s2 
of  the  samples  cut  off  in  the  direction  orthogonal  to  the  maximal  defonnation  processed  in 
Experiment  #5. 

Potential  for  “net  shape”  HIP  technology  based  on  advanced  process  modeling 
accounting  inherent  process  anisotropy,  to  reduce  cost  and/or  time  required  to  produce  specific 
liquid  rocket  engine  components  was  analyzed. 

This  development  effort  has  provided  a  general  proof  of  viability  of  proposed 
“net  shape”  HIP  technology  based  on  the  advanced  process  modeling. 
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